Phase Diagrams of the Two-Orbital Hubbard Model with Different Bandwidths 
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We investigate the two-orbital Hubbard model in infinite dimensions by means of the self-energy 
functional method. By calculating the entropy, susceptibility, and quasi-particle weight at zero 
temperature, we determine the phase diagram for the system with same and difi'erent bandwidths, 
which is compared with that obtained recently. It is clarified that orbital fluctuations play a key 
role in controlling the nature of the Mott transitions in the system. 
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I. INTRODUCTION 



Transition metal oxides have attracted much interest 
in their various aspectspi^^^ Among them, heavy fermion 
behavior in d-electron systems stimulates experimen- 
tal and theoretical investigations. One of typical ex- 
amples is the lithium vanadate LiV204^ where heavy 
fermion behavior with the enhanced specific-heat coeffi- 
cient 7e — 0A2J/molK^ has been observed at low tem- 
peratures. It has been suggested that the large mass 
enhancement in this system may originate from geomet- 
rical frustration. "^-^-^ Recently, it has been pointed out 
that degenerate orbitals in the t2g subshell also play an 
important role in realizing heavy fermion behavior in 
LiV204.^'^ Another example is the isovalent ruthenate 
arroy Ca2_xSrxR'U04. In the compound, the substitu- 
tion of Ca^+ ions for x > 0.5 realizes heavy fermions 
in the t2g subshell,^ where some of physical quanti- 
ties have a cusp singularity around x — 0.5. Further- 
more, the unexpected s — 1/2 moment per Ru-ion co- 
existing with the metallic state has been observed in 
that region,^ which suggests the existence of the orbital- 
selective Mott transition (OSMT)iSiii This type of the 
Mott transition has also been suggested to occur in the 
compound Lan-nNinOsn-n /^'^'^ stimulating further ex- 
perimental investigations on the multiorbital systems. 

In the above compounds, electron correlations 
in the system with orbital degeneracy are impor- 
tant in understanding heavy fermion behavior. In 
particular, the concept of the OSMT provides a 
new paradigm of metal-insulator transitions in the 
multiorbital systems. Although intensive theo- 
retical studies on these topics have recently been 



with different bandwidths. It has been pointed out 
that the double OSMTs occur in general, which merge 
to a single Mott transition only under the special 
condition. On the other hand, it has been 
proposed that the double transitions always occur if the 
system has a large difference of the bandwidthsi^2i2S 
Therefore, it is highly desired to discuss the nature of 
the Mott transition in order to clarify what triggers the 
single or double transitions in the system with different 
bandwidths. 

Motivated by this, we investigate the Mott transitions 
in the multiorbital systems. In particular, we focus on 
orbital fluctuations in this paper, which may play an 
important role in stabilizing the metallic state in the 



system 
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By making use of the self-energy functional 



approach (SFA) proposed by PotthofffSiS we confirm 
some previous works^^i^ and determine the phase di- 
agrams. We examine how orbital fluctuations are af- 
fected by the Hund coupling and/or the difference of the 
bandwidths to clarify that the enhanced orbital fluctua- 
tions make single Mott transition stable against double 
OSMTs. 

This paper is organized as follows. In Sec. m we in- 
troduce the two-orbital Hubbard model and briefly sum- 
marize the SFAiSiS The nature of the Mott transitions 
in the two-orbital system with same bandwidths is dis- 
cussed in Sec. IIIII We also discuss how the OSMT is 
realized in the system with different bandwidths in Sec. 
IIVI A brief summary is given in Sec. 



there are some controversial conclusions even in the 
simplest multiorbital model. In fact, in the infinite- 
dimensional two-band system with equivalent orbitals, 
the existence of the Hund coupling induces the first- 
order Mott transition to the insulating phase>ii22ii4 
However, recent results obtained by dynamical mean- 
field theory (DMFT)^°i^^ combined with the numerical 
renormalization group claim that the Mott transition 
is of second-order in a certain parameter region^ In 
addition, there are some open questions in the system 



II. MODEL AND METHOD 
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We consider the two-orbital Hubbard model with dif- 
ferent bandwidths, which is given by the following Hamil- 
tonian, H ^ Ho + W , W = K with 



^0 = 51 --"%)cLaCia<T, (1) 

<i,j> ,a,cr 

-J(c]'^Cjii42|Ci2T + c]'^c|^C^2TC^21 + H.C.), (2) 
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where cl^^{ciaa) creates (annihilates) an electron with 
spin cr(=t: i) E^iid orbital a{— 1,2) at the i th site, and 
Uiaa is the number operator. Here, ta denotes the hop- 
ping integral for the ath orbital, /i the chemical potential, 
U{U') the intra-orbital (intcr-orbital) Coulomb interac- 
tion, and J the Hund coupling including the spin-flip and 
pair-hopping terms. In the paper, we impose the condi- 
tion U — U' + 2J due to the rotational symmetry of 
degenerate orbitals. 

To discuss the Mott transitions in the multior- 
bital systems, we make use of the SFA. Since this 
method is based on the variational principle, it has 
an advantage to discuss the properties of the Mott 
transition systematically. In fact, it has success- 
fully been applied to various systems such as Hub- 
bard models in infinite dimensions^SiS^iS^iSi^ and finite 

dimensionS|57i58,59i60,61,62,63 

In this approach, the ground potential Q is given as. 



nm = F[S] + Trln[-(G, 
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(3) 



where PI'S] is the Legendre transformation of the 
Luttinger-Ward potential.^* Gq and S are the bare 
Green function and the self-energy, respectively. Here 
the Dyson equation G^^ = Gg"^ — S is obtained un- 
der the condition ari[S]/aS = 0^ where G is the full 
Green function. We wish to note that the potential F[E] 
does not depend on the detail of the non-interacting 
HamiltonianiS This allows us to introduce a reference 
system with the same interacting term. The Hamilto- 
nian is exphcitly given by Href (f) = 'Wo(t') -I- H' with 
the parameter matrix t'. Then we obtain the grand po- 
tential as, 

f][s(t')] = n(t') 

+ Trln[-(cj + /i-t-S(t'))"^] 

- Trln[-(c. + /i-t'-S(t'))-i], (4) 

where f2(t') and S(t') are the grand potential and the 
self-energy for the reference system. The condition 
dQ['E{t')]/dt' = gives us an appropriate reference sys- 
tem Tiref (f) in the framework of the SFA. 

If the SFA is applied to infinite-dimensional correlated 
electron systems, an Anderson impurity model is one of 
the most appropriate reference systems, which is given 
by 



nAi) 



(5) 



Nt, 



(i)„Wt„(i) 



fe=l aa 



EE^(4.«i2.+^-)+w:, (6) 



k=l aa- 



where O'ki^i'^kaa) creates (annihilates) an electron with 
spin (7 and orbital a at the k{= 1, 2, • ■ • -/Vb)th site, which 



is connected to the ith site in the original lattice. In 
the limit of Nt oo, the condition 9f][S(t')]/9t' = 
reproduces self-consistent equations of DMFTi^SiSi^ 

The grand potential per site is rewritten as, 



n/L = a.„p-2^^F(^:,J-2^5]i^(4j 



am a. k—1 

n — QO 

+ 2^^ / dzp^{z)F[LUa,m{z)], (7) 

F{x) = -Tln[l + exp(-a;/T)] (8) 

where ilimp is the grand potential for the reference sys- 
tem. ujam{z) [ijj'am] ^hc polc of Green function G {G') 

£ka — 



for the original (reference) system and 
The Green function of the reference system is given as, 

G'^iuj) = [c^ + /i-eoa- A„(w)-E„(w)]^i, (9) 



kOL 



(10) 



fe=l 



'ka 



where Sq(w) is the self-energy for the ath orbital. On 
the other hand, the full Green function is given as. 



(11) 



Note that the following differential equation is efficient 
to deduce the poles ujamiz) in the Green function. 



dz 



(12) 



By solving these equations, we estimate the grand poten- 
tial numerically to discuss the effect of electron correla- 
tions. 

To clarify the nature of Mott transitions, we calculate 
various physical quantities. In the metallic phase, the 
quasi-particle weight for the ath orbital, Za^ is useful to 
discuss how the Fermi liquid states are renormalized by 
the Coulomb interactions. This quantity is proportional 
to the inverse of the effective mass, which is defined as. 



= 1- 



du! 



(13) 



Furthermore, we calculate the local susceptibilities, 
which may clarify how spin, orbital and charge fluctu- 
ations affect the stability of the metallic state. These 
read 



Xs = 



Xo 



(rK(r)-nx(r)][nT(0)-nx(0)])dr, 



(r[ni(T) - n2(r)][ni(0) - n2(0)])dr. 



{T[n{T)-2][n{Q)-2])dT, (14) 
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where n = J^aa'^aa, n^{a) = Y.a{a)'^aa and 

naa = cl^^Ciaa, ^ the time ordered operator, A{t) = 
q-'Htj^^Ht g^j^jj p |.jjg inverse teraperature. Note that in 
infinite dimensions, these local quantities coincide with 
those for the reference system with iVf, — > ooi^ In the pa- 
per, calculating these quantities approximately in terms 
of the reference system with a finite Ni,, we discuss how 
fluctuations affect the stability of the metallic state. We 
also calculate the residual entropy in the system to char- 
acterize the Mott insulating phase, where localized elec- 
trons are realized with the free spins. It is given as. 



S = - — 



dT 



(15) 



= S; 



imp 



a k=l 



2EE 



dzpa{z)F'[uJarn{z)], 



F'{x) 

fix) 



[Fix) ~{x- dx/dT)fix)]/T, 
l/(l + exp(x/r)). 



(16) 

(17) 
(18) 



where S'imp = —d^hnp/dT is the entropy of the refer- 
ence system and we have used the condition dVl/dt' = 0. 
Note that the temperature derivative of the poles for the 
reference system is zero: duj'^^/dT = 0. By estimating 
dojamiz) / dT carefully,— we can calculate the entropy for 
the system. 

Here, we use the semicircular density of states Pa(i^) — 
A/nWa^l - {2x/Wa)^, where Wai= 4i„) is the band- 
width for the ath orbital, which corresponds to an infinite 
coordination Bethe lattice. In this paper, we restrict our 
discussions to the paramagnetic case in the half-filled sys- 
tem, by setting the chemical potential /x = U/2+U' — J/2. 
In the following, by varying the ratio of the bandwidths 
R = Wi/W2{< 1) with a fixed ^2 = 1 (that is an energy 
unit), we proceed to discuss the Mott transitions in the 
two-orbital model. 



III. MOTT TRANSITION IN THE SYSTEM 
WITH SAME BANDWIDTHS 

Let us consider the two-orbital system with same 

bandwidths. Mott transitions in the system have been 

discussed so far, by combining DMFT with numerical 
techniques 10-ii-i4.15.i6.17.i8.19.20.21.22.23.24.25.26.27.28.29.30.31 



We first calculate the ground potential as a function 
of = Vi 1 = Vi 2 at zero temperature, as shown in Fig. 

12 



0.04 



o 

a 



-0.04 



-0.08 



o 

a 



0.04 1 


-0.04 
-0.08 



j/u=o ^^^^^ 


7/{7=0.01 


u=% 

i7=8.5 

— U=9.1 
--- {7=9.5 


(7=6 \ 

(7=6.5 

. — (7=6.91 
— ■ (7=7.5 


\ \ r- 

//C/=0.03 


\ \ 7 — 

//{7=0.25 


" {7=5 \ 

{7=5.5 \ 

— {7=5.84 


(7=3 " 

(7=3.2 

— (7=3.52 

.... [;^4 



0.2 



0.4 



0.2 



0.4 



An exact diagonalization study claimed that the first- 
order transition occurs with the jump in physical 
quantitiesi^^ In contrast, it was suggested that the 
second-order transition occurs in a certain parameter re- 
gion by means of the numerical renormalization group,-^? 
In this section, to resolve the controversial conclusions 
for the properties of the Mott transition, we make use 
of the SPA with A^f, = 1. At the end of the section, 
we check the validity of our analysis by comparing the 
results of A^b = 1 with those of A;, = 3. 



FIG. 1: The grand potential as a function of the variational 
parameter V . Thick lines correspond to the critical value (7c. 

When U is small, it is found that the grand potential 
has a minimum at finite V . The hybridization Y in the 
reference system, roughly speaking, represents the effec- 
tive bandwidth for the original system. Therefore, the 
metallic state is stabilized in this case. It is seen that 
the increase of the Coulomb interaction with a fixed ra- 
tio 3 jU = shifts the stationary point toward the origin 
continuously. This implies that the effective bandwidth 
is gradually decreased, and the transition occurs to the 
Mott insulating phase at C/ = J7c- In contrast, the intro- 
duction of the Hund coupling leads to different behavior. 
For instance, we focus on the system with the weak Hund 
coupling J jU — 0.03. When J7 = 5, a minimum appears 
at y 0.58, where the metallic ground state is realized. 
The increase of the Coulomb interaction induces another 
minimum around ~ 0, which represents the insulating 
state. This double-well structure in the grand potential 
suggests the existence of the first-order Mott transition in 

,32.S^34„^5,SR.^.3K3B.AB.^j4a^3.A4A5^6.47;..lg,49F^. gach State, WC 

determine the critical value Vc = 5.84. A similar struc- 
ture is also observed in the case J/U — 0.25. However, 
the potential barrier between these two states becomes 
small at the first-order transition point, as shown in Fig. 
n Therefore, it is expected that the singularity charac- 
teristic of the first-order transition is more difficult to be 
observed in this case. 

These features mentioned above can be seen in the 
physical quantities. Here, we show the entropy per site 
S/L in Fig. |21 When J/U = 0, the increase of the 
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FIG. 2; The entropy per site as a function of U . 



Coulomb interaction triggers the Mott transition aXU — 
Uc, where the entropy jumps up to S/L = In 6. The 
value is explained by the fact that sixth degenerate states 
are lowest in the atomic limit. On the other hand, the 
Hund coupling forms the triply degenerate states in the 
limit. Therefore, when J/U ^ 0, S/L = ln3 appears in 
the insulating phase. It is found that the second-order 
transition occurs in the system J/U = 0, while the first- 
order transition with the hysteresis occurs in the system 
with the Hund coupling. For example, when J/U = 0.03, 
the metallic state is stabilized up to the critical point 
t/*^ — 6.12, while the insulating state is stabilized down 
to the critical point U/. = 4.82. Note that these critical 
points J7*^ and U/. induced by the Hund coupling at zero 
temperature are different from critical points Ud and 
Uc2^ which are defined as Mott critical points at finite 
temperatures. Therefore, an inequality U/. < U//^ is not 
necessarily satisfied, which will be discussed in detail in 
the next section. 

By performing similar calculations, we obtain the 
phase diagram as shown in Fig. O We find that the 
phase boundary C/*^ is always larger than U/. except for 
the special condition J/U = 0. This implies that the 
first-order Mott transition, in general, occurs in the two- 
orbital system with same bandwidths. When the sys- 
tem has the special condition J/U = 0, the critical point 
U/. (= U^ = 9.24) is different from the extrapolation of 
the curve U/. 44 This exception may originate from the 
discontinuity of the residual entropy for the Mott insu- 
lating phase. It is also found that the coexisting region 
bounded by these two lines shrinks with the increase of 
the Hund coupling J. This interesting feature may be 
due to orbital fluctuations, which sometimes play an im- 
portant role in understanding the Mott transitions in the 
two-orbital system. 

To make this point clear, we also calculate the local 
orbital susceptibility in the metallic state [U < U^'^), as 
shown in Fig. ^ When J/U = 0, the increase of the 
Coulomb interaction enhances the spin and orbital sus- 



FIG. 3: The zero-temperature phase diagram for the degener- 
ate Hubbard model with same bandwidths R = 1. U' iU^') 
is the critical point where the insulating (metallic) state dis- 
appears. Uc represents the first order transition points and 
the filled region represents the coexisting state. 
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FIG. 4: The local orbital susceptibility in the system with 
same bandwidths. The inset shows orbital and charge sus- 
ceptibilities in the case J = 0.0317. 

ceptibilities (the spin susceptibility is not shown). Even- 
tually, these quantities diverge simultaneously at the 
transition point. On the other hand, when one turns 
on the Hund coupling [J/U = 0.01), different behavior 
appears. The introduction of the Coulomb interaction 
first enhances orbital fluctuations, by reflecting the high 
symmetry U ^ U' . Further increase of the interactions 
renormalizes electrons in both orbitals, where the Hund 
coupling has a tendency to suppress orbital fluctuations. 
In fact, in contrast to the case U = U', the orbital suscep- 
tibility is much smaller around the Mott critical point. 
We also find the sudden decrease of the orbital suscepti- 
bility where the first-order Mott transition occurs. When 
the system is located far from the condition [/ = [/', the 
increase of the Coulomb interaction suppresses orbital 
fluctuations strongly. Therefore, a tiny jump appears in 
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the curve at = 3.54 {J/U = 0.25). 

The obtained results may shed hght on the nature of 
the Mott transitions in the two-orbital systems with same 
bandwidths. When J/U = 0, the second-order Mott 
transition occurs to the insulating phase. On the other 
hand, the introduction of the Hund coupling induces the 
first-order transition with the hysteresis in the physical 
quantities. It is also found that as J is increased, orbital 
fluctuations are suppressed strongly, where the singular- 
ity for the transition becomes obscure. In this case, the 
system can be regarded as the system with independent 
orbitals. Therefore, the OSMT is expected to occur if 
the bandwidths are different from each other, which will 
be discussed in the next section. We wish to comment 
on the related recent work.^^^ The introduction of the 
Hund coupling is known to induce the first-order tran- 
sition, while the nature of transitions in the large J re- 
gion has not been clarified up to now. In the region, 
the first- (second-) order transition was suggested by the 
exact diagonalization (numerical renormalization group) . 
On the other hand, our systematic analysis shows the 
jump singularity in the curve for each physical quantity, 
except for the J — case. Therefore, we believe that fur- 
ther increase of the Hund coupling does not induce the 
second-order Mott transition in the system with same 
bandwidths. 

Before closing this section, we would like to check the 
validity of our analysis. In Fig. (SJ we show the quasi- 
particle weights obtained by the SFA with Ni, = 1 and 
Nf, = 3. If the Coulomb interactions are introduced. 
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U 

FIG. 5: The quasiparticle weight Z as a function of U ob- 
tained from the reference systems Nh = 1 and A*';, = 3 when 
J/U = 0. 

the quasi-particle weight is decreased from unity. In 
this region, it is found that the quasi-particle weight ob- 
tained for the system A^;, = 1 is slightly larger than that 
for the system Nt, — 3. On the other hand, when the 
system approaches the Mott transition point, the quasi- 
particle weight is hardly affected by the number of the 
sites for the reference system. This result may be ex- 



plained by the fact that the simplified reference system 
{Nf, — 1) can describe the low-energy physics around 
the Fermi level and thereby we can determine the crit- 
ical point quantitativelyi^SiSi^ However, the reference 
system A^;, = 1 does not describe the high energy part 
properly, resulting in the overestimate of Z in the inter- 
mediate region. 

Similar behavior is also observed in the finite J case, as 
shown in Fig. El As the Hund coupling is increased, the 




FIG. 6: The quasiparticle weight Z as & function of U in the 
cases J/U = 0.01 and J/U = 0.25. 

Fermi liquid state is stabilized up to J7*^ ~ 7.3 and 3.5 at 
J/U = 0.01 and 0.25 respectively, where the first-order 
transition occurs. It is found that the jump singularity 
for Aft = 3 is weaker than that for A^;, = 1, which may 
expect that the transition becomes of second order in the 
Aft — > oo case. According to the results for the single- 
band model)SSiSi^ however, the physical quantities for 
A^b = 3 nearly equal those for A'b > 3. Therefore, we 
believe that the nature of the Mott transition is hardly 
affected by A^, and the first-order transition remains in 
the A^f, oo case. In addition, we also find that the 
critical point U^^ weakly depends on the number of sites 
Nb, implying that the obtained phase boundaries shown 
in Fig. 13 have been determined rather well even by the 
reference system Nf, = 1. 



IV. ORBITAL-SELECTIVE MOTT 
TRANSITIONS VS. SINGLE MOTT 
TRANSITION 

In this section, we consider Mott transitions in the 
two-orbital system with different bandwidths (i? ^ 

j^-j 10.11.31.32.33.34.35.36.37.38.39.40.41.42.43.45.46.47.48.49 ^las 

been suggested that in the system with R = 0.5, the 
OSMT occurs in the case J/U = 0.25, while a single 
Mott transition occurs in the case J/U = Q^^iSk On the 
other hand, it was recently claimed that the double tran- 
sitions always occur in the small R case M/—- By means 
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of the SFA with Nf, = 1, we determine the detailed phase 
diagram to discuss the nature of the Mott transitions 
systematically. 

First we consider the system without Hund coupling 
{J/U = 0). In the system with same bandwidths {R = 
1), the second-order Mott transition occurs at Uc — 9.24, 
where the quasi-particle weight continuously reaches zero 
and the orbital susceptibility diverges. Away from the 




FIG. 7: The quasi-particle weights for a th orbital Za when 
ii = 1 and 0.2. 




FIG. 8: The inverse of orbital and spin local-susceptibilities 
as a function of U when J/U = 0. The inset shows orbital 
and charge susceptibilities, and the entropy per site when 
R^O.17. 

condition (e.g. R — 0.2), the introduction of the in- 
teraction decreases the quasi-particle weights Zi and Z2 
in different ways reflecting the difference of bandwidths, 
as shown in Fig [7| This enhances the spin fluctuations 
and suppresses orbital fluctuations. However, when the 
system approaches the Mott transition point, the effect 
of the bare bandwidths is diminished due to the strong 
renormalization, which enhances both spin and orbital 
fluctuations. Therefore, in the case R = 0.2, monotonic 



(non-monotonic) behavior appears in the curves of the 
spin (orbital) susceptibility. As a consequence, a single 
Mott transition occurs, where the quasi-particle weight 
for each orbital vanishes simultaneously.'^^ On the other 
hand, somewhat different behavior appears in the case 
R < Rc{^ 0.192), as shown in the inset of Fig. El If the 
Coulomb interactions are increased, orbital fluctuations 
are strongly suppressed due to the difference of the effec- 
tive Coulomb interactions. In this case, the orbital sus- 
ceptibility never diverges, which suggests the existence of 
the OSMT, as discussed by Medici et alm^ and Ferrero et 
alM: In fact, the OSMT yields a locahzed spin s = 1/2 at 
each site. Therefore, the residual entropy S/L — ln2 ap- 
pears and the orbital susceptibility merges to the charge 
susceptibility at ^ 3.2, as shown in the inset of Fig. 
|S1 Further increase of the Coulomb interaction decreases 
the charge and orbital susceptibilities and another Mott 
transition occurs at ~ 5.9. 

By estimating the critical points systematically, we ob- 
tain the phase diagram shown in Fig. O When R = 1, 
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FIG. 9: The phase diagram for the case J/U = 0. U' {U^') 
is the critical point where the insulating (metallic) state be- 
comes unstable with increasing (decreasing) U. The dotted 
lines denote the phase boundaries obtained by the simple es- 
timation (see text). 

the model is reduced to the system with same band- 
widths, where the Mott transition occurs between the 
metallic and Mott insulating phases. Even when the sys- 
tem has two distinct orbitals, the single Mott transition 
still occurs {U^^ = U^) owing to enhanced orbital fluc- 
tuations within the region {R^ < i? < 1). In this region, 
the phase boundary is simply deduced by the Gutzwiller 
approximation as, Uc oc J deep{e) = (1 -I- i?)/2, which 
is in good agreement with the results obtained by the 
SFA. On the other hand, in the region R < R'^, the dou- 
ble Mott transitions occur (J7*^ < U^) when the inter- 
action is varied. Then the orbital-selective Mott insu- 
lating (OSMI) phase appears between the metallic and 
insulating phasesi^SiSS In the case, the system may be 
regarded as two independent single-band systems. In 
fact, the phase boundaries ~ UcR and U^' ^ Uc, 
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where Uc = 5.85 is the critical point for the single-band 
model^^^iS^ are in good agreement with our SFA results 
in the case i? ^ 1, as shown in Fig. The phase dia- 
gram is consistent with that obtained by Medici et alm^ 
and Ferrero et aL^^. 

Next, we consider the effect of the Hund coupling in 
the system with different bandwidths R = 0.5. In Fig. 
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FIG. 10: The local susceptibilities and the entropy per site in 
the system with different bandwidths R = 0.5. 

1101 we show the physical quantities as a function of U 
with a fixed ratio J/U. When the Hund coupling is 
small J/U = 0.01, the interaction rcnormalizes electrons 
depending on each bandwidth, and it triggers the first- 
order single transition to the Mott insulating phase at 
U^^ = 5.47. Then, the singularity appears in the curve 
of the susceptibility and the entropy. This implies that 
the existence of the Hund coupling induces the first-order 
Mott transition with the hysteresis, where the Mott in- 
sulating phase is stabihzed down to C// = 4.72. Al- 
though this result is the same as that in the system with 
i? = 1, it does not necessarily imply that the difference 
of the bandwidths is irrelevant in the whole parameter 
space. Indeed, further increase of the Hund coupling sup- 
presses orbital fluctuations, leading to double transitions 
in the system. When J/U = 0.1, the introduction of 
the Coulomb interaction decreases the orbital and charge 
susceptibilities. Consequently, these quantities merge to 
each other and the spin susceptibility diverges at the crit- 
ical point U^'^ = 3.34. This means that the second-order 
OSMT occurs in the narrower band. Thus, free spins 
s — 1/2 are induced at each site, yielding the residual 
entropy S/L = In 2. Furthermore, as the interactions are 
increased, both orbital and charge susceptibilities reach 
zero at U^ = 4.06. 

Performing similar calculations with several choices of 
R, we end up with the phase diagrams shown in Fig. 1111 
An important point is that the phase boundaries U^'^ and 
[// cross each other at i? = Rc, where Rc = 0.30 and 0.82 
for J/U = 0.1 and 0.01, respectively. When R < Rc, the 
difference of the bare bandwidths is essential to stabilize 
the OSMI phase between two phase boundaries {U^^ < 
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FIG. 11: The phase diagram when J/U = 0.01 and J/U = 
0.1. The gray region denotes the coexisting state. 



U < Ul). However, when R > Rc, the metallic phase 
and the Mott insulating phase coexist in the region U/. < 
U < Uc'^ . In this case, the nature of Mott transitions 
is similar to that in the system with same bandwidths, 
where the first-order Mott transition occurs. The above 
result may suggest that the existence of the OSMT is 
associated with the first-order transition in the multi- 
orbital Hubbard model with same bandwidths. 

To discuss how the coexisting phase competes with the 
OSMI phase, we also show the phase diagrams with a 
fixed R in Fig. [HI When R < 0.192, there always 




FIG. 12: The phase diagram for the degenerate Hubbard 
model with different bandwidths R = 0.1, 0.5 and 0.8. The 
black and gray region denote the OSMI phase and the coex- 
isting state, respectively 

appears the OSMI phase due to the large difference of 
bandwidths. Note that the the phase boundary some- 
times has a discontinuity at J/U = 0, which may orig- 
inate from the spin and orbital degeneracy of the Mott 
insulating ground state, as discussed in the previous sec- 
tion. In contrast, the OSMI phase and the coexisting re- 
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gion appear in the case 0.192 < i? < 1. When the system 
has the large Hund couphng, orbital fluctuations are sup- 
pressed and the OSMI phase is stabilized. The decrease 
of J enhances orbital fluctuations in the metallic phase, 
making the OSMI phase unstable. The phase boundaries 
then merge at the critical point. Further decrease of the 
Hund coupling induces a first-order single Mott transi- 
tion due to enhanced orbital fluctuations, where the co- 
existing phase appears instead of the OSMI phase. When 
i? — > 1, the difference of the bandwidths becomes irrel- 
evant. Therefore, the region of the OSMI phase shrinks 
and the coexisting phase becomes stable. 

Up to now, there were some controversial conclusions 
for the nature of the OSMT. Most of studies claimed that 
the second order OSMT occurs in generali^Siii In con- 
trast, Medici et al.^^ and Arita et alm^ recently claimed 
that the OSMT is of first order by combining DMFT 
with the exact diagonalization and the projective quan- 
tum Monte Carlo simulations. In our paper, we have 
calculated the grand potential in the system carefully to 
clarify that the OSMT is of second order. Although the 
accuracy of our calculations depends on the size of the 
reference system in the framework of the SFA, the na- 
ture of Mott transition is hardly affected by the size, as 
discussed before. Therefore, we believe that the obtained 
results may shed light on the nature of the Mott transi- 
tions in the two-orbital systems. 

V. SUMMARY 

We have investigated the two-orbital Hubbard model 
in infinite dimensions to discuss the nature of Mott tran- 



sitions at half-filling. By making use of the self-energy 
functional approach, we have clarified that the introduc- 
tion of the Hund coupling induces the first-order tran- 
sition in the system with same bandwidths. It has also 
been found that the increase of the Hund coupling sup- 
presses orbital fiuctuations, where the singularity char- 
acteristic of the first-order transition becomes obscure. 
Namely, the upper and lower critical points for the tran- 
sition approach each other closely, where the weak first- 
order transition occurs. This behavior is also affected 
by the difference of the bandwidths. In fact, the differ- 
ence of the bandwidths suppresses orbital fluctuations, 
making the first-order transition unstable. As a conse- 
quence, the double second-order transitions occur, where 
the OSMI phase appears between the metallic and the 
Mott insulating phase. Our obtained results reproduce 
some previous works on the system with same and dif- 
ferent bandwidths. ^SiSSiSSi^ Furthermore, taking into ac- 
count of orbital fluctuations carefully, we have resolved 
some controversial conclusions for the nature of Mott 
transitions at half-filling. 
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